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Abstract. We introduce concepts of "recursive polynomial remainder sequence (PRS)" and 
"recursive subresultant," and investigate their properties. In calculating PRS, if there exists 
the GCD (greatest common divisor) of initial polynomials, we calculate "recursively" with 
new PRS for the GCD and its derivative, until a constant is derived. We call such a PRS a 
recursive PRS. We define recursive subresultants to be determinants representing the coef- 
ficients in recursive PRS by coefficients of initial polynomials. Finally, we discuss usage of 
recursive subresultants in approximate algebraic computation, which motivates the present 
work. 



1 Introduction 

The polynomial remainder sequence (PRS) is one of fundamental tools in computer algebra. Al- 
though the Euclidean algorithm (see Knuth ([!]) for example) for calculating PRS is simple, co- 
efficient growth in PRS makes the Euclidean algorithm often very inefficient. To overcome this 
problem, the mechanism of coefficient growth has been extensively studied through the theory of 
subresultants; see Collins ([2]), Brown and Traub ([3j), Loos ([!]), etc. By the theory of subresul- 
tant, we can remove extraneous factors of the elements of PRS systematically. 

In this paper, we consider a variation of the subresultant. When we calculate PRS for poly- 
nomials which have a nontrivial GCD, we usually stop the calculation with the GCD. However, 
it is sometimes useful to continue the calculation by calculating the PRS for the GCD and its 
derivative; this is necessary for calculating the number of real zeros including their multiplicities. 
We call such a PRS a "recursive PRS." 

Although the theory of subresultants has been developed widely, the corresponding theory for 
recursive PRS is still unknown within the author's knowledge: this is the main problem which we 
investigate in this paper. By "recursive subresultants," we denote determinants which represent 
elements of recursive PRS by the coefficients of initial polynomials. 

This paper is organized as follows. In Sect. [2 we introduce the concept of recursive PRS. In 
Sect.[3j we define recursive subresultant and show its relationship to recursive PRS. In Sect. [4] we 
discuss briefly using recursive subresultants in approximate algebraic computation. 



2 Recursive Polynomial Remainder Sequence (PRS) 

First, we review the PRS, then define the recursive PRS. At last, we show recursive Sturm sequence 
as an example of recursive PRS. 



2.1 Definition of Recursive PRS 

Let R be an integral domain and polynomials F and G be in R[x] . We define a polynomial remainder 
sequence as follows. 

* This research was partially supported by Japanese Ministry of Education, Culture, Sports, Science and 
Technology under Grant-in-Aid for Young Scientists (B), 14780181, 2002. 
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Definition 1 (Polynomial Remainder Sequence (PRS)). Let F and G be polynomials in 
R[x] of degree m and n (m > n), respectively. A sequence 

(Pu---,Pi) (1) 

of nonzero polynomials is called a polynomial remainder sequence (PRS) for F and G, abbreviated 
to prs(P, G), if it satisfies 

Pi =F, P 2 = G, OiPi-2 = qt-iPi-i + Pi Pi, (2) 

for i = 3, . . . , I, where a^, . . . , a>i, 0s, . . . ,0i are elements of R and deg(Pi_i) > deg(Pj). A sequence 
((03, 0s), ■ ■ ■ ! (o-hPi)) is called a division rule for prs(F, G) (see von zur Gathen and Lucking (®))- 
If Pi is a constant, then the PRS is called complete. □ 

If F and G are coprime, the last element in the complete PRS for F and G is a constant. 
Otherwise, it equals the GCD of F and G up to a constant: we have prs(F, G) = (Pi = F, P2 = 
G, . . . ,Pi = 7 • gcd(P, G)) for some 7 G R. Then, we can calculate new PRS, prs(P u £Pi), and if 
this PRS ends with a non-constant polynomial, then calculate another PRS for the last element, 
and so on. By repeating this calculation, we can calculate several PRSs "recursively" such that the 
last polynomial in the last sequence is a constant. Thus, we define "recursive PRS" as follows. 

Definition 2 (Recursive PRS). Let F and G be the same as in Definitional^ Then, a sequence 

fpU) p(l) p(2) p(2) p(t) pit)-, /ox 

l/l ' * * ' ' Zl ' 1 ' ■ ■ • I 'la 1 • ■ ■ 1 r \ ' • • • ! r i, I \°) 

of nonzero polynomials is called a recursive polynomial remainder sequence (recursive PRS) for F 
and G, abbreviated to rprs(i 7 ', G), if it satisfies 

P[ 1] =F, P 2 (1) =G, Pi 1) =7i-gcd(P 1 (1) ,P 2 (1) ) with 7l eR, 
( pW,P«,...,P«) = prs(P«,P«), 

Pf^Pt 1 ', H k) = -£Pt?> P^=l k -Scd(Pt\P 2 {k) ) wtih lk tR, ^ 
(pf\P 2 (&) ,...,p/ fc fe) )=prs(P 1 (fc) ,P 2 (fc) ), 
fork = 2,...,t.If af ] , # } € P satisfy 



a 



/or fc = 1, . . . , t and i = 3, ...,lk, then a sequence {(cx^ , 0^), ■ ■ ■ , (a|*\ /?|*')) is called a division 
rule for rprs(P, G). Furthermore, «/P; is a constant, then the recursive PRS is called complete. 



□ 



2.2 Example of Recursive PRS: Recursive Sturm Sequence 

Sturm sequence is a variant of PRS, which is used in Sturm's method, for calculating the number 
of real zeros of univariate polynomial (for detail, see Cohen ([7]) for example). Note that Sturm's 
theorem is valid for not only polynomials having simple zeros but also those having multiple zeros 
(see Bochnak, Coste and Roy ([5]) for example). Here, we define "recursive Sturm sequence" to 
calculate the number of real zeros including multiplicities, as follows. 

Definition 3 (Recursive Sturm Sequence). Let P(x) be a real polynomial of degree m. Let a 
sequence of nonzero polynomials be defined by a recursive PRS in Definition^ calculated as 

(complete) rprs(P(x), — P(x)), (6) 

with division rule given by 

(af\/f>) = (l,-l), (7) 

for k — 1, . . . , t and i = 3, . . . , Ik- Then, the sequence (J6j is called the recursive Sturm sequence of 
P{x). □ 
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Example 1 (Recursive Sturm Sequence). Let P(x) = (x + 2) 2 {(x — 3)(x + l)} 3 , and calculate 
the recursive Sturm sequence of P(x). The first sequence Li = (P^\ P 4 (1) ) has the following 
elements: 

P« = P( x ) = ( x + 2) 2 {(x - 3)(x + l)} 3 , 

= ^-P(x) = Sx 7 - Ux 6 - 102a; 5 + 80a; 4 + 460a; 3 + 66a; 2 - 558a; - 324, 
ax 

/i) 75 e 45 = „„ 4 225 , 3315 2 4815 945 ( 8 ) 
„m 128 * 256 4 256 , 1024 2 4224 2304 

p * = ^ - 25^ - ~ x + ^r x + ^r x+ ^r- 

(2) (2) 

The second sequence L 2 = (P{ , . . . , P 4 ) has the following elements: 

p (2) _ p (i) _ 128 8 256 4 256 3 1024 4224 2304 
A _ ~25" X ""25* + ^5~ a: + ^5~ X+ ^5~' 

d „m 128 „ 1024 , 768 2048 4224 



da; 4 5 25 5 25 25 

(2) _ 14848 3 1536 2 88576 66048 

3 625 125 625 625 
(2 ) _ 12800 2 25600^ 38400 

4 841 X 841 X ~ 841 

The last sequence L3 = (Pf 3 \ . . . , P^) has the following elements: 

,(3) _ p (2) _ 12800 2 25600 38400 



(9) 



841 841 841 ' 

w - w = 25600^25™ (1Q) 



,0) _ _d (2) _ 25600 _ 25600 
2 ~ dx 4 ~ 841 X 841 ' 
( 3) _ 51200 



p 

3 841 

For PRS ife, k = 1,2, 3, define sequences of nonzero real numbers A(Lfc, — exo) and A(Lfc, +00) as 

A(L fc ,-oo) = ((-lK , lc(P 1 W ),...,(-l)< ) lc(pf )), 
A(i fc ,+oo) = (lc(P 1 (fe) ),...,lc(p( fc) )), 

where = deg(P^) denotes the degree of P/ fe ' ) and lc(P^) denotes the leading coefficients of 
Then, A(ifc, — 00) and A(Lfc,+oo) for fc = 1,2,3 are 

75 128 
A(Li,±oo) = (l,±8,-,±— ), 

y ' ' K 25 ' 5 ' 625 ' 841 ; ' ' 
, ,12800 25600 51200, 

For a sequence of nonzero real numbers L = (ai, . . . , a m ), let V^L) be the number of sign variations 
of the elements of L. Then, we calculate the number of the real zeros of P(x), including multiplicity, 
as 

3 

Y,{V{X{L k , -00)) - V(\(L k , +00))} = 3 + 3 + 2-8. (13) 

□ 
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3 Subresultants for Recursive PRS 

Let F and G be polynomials in R[x] such that 

F(x) = f m x m + ■■■ + fox , G(x) = g n x n + ■ ■ ■ + g a x°, with m > n > 0. (14) 

To make this paper self-contained and to use notations for our definitions, we first review the 
fundamental theorem of subresultants, then discuss subresultants for recursive PRS. 



3.1 Fundamental Theorem of Subresultants 

There exist several different definitions of subresultants. Here, we follow definitions and notations 
basically by von zur Gathen and Lucking ([6]). 

Definition 4 (Sylvester Matrix). Let F and G be as in (jT4j) . The Sylvester matrix of F and 
G, denoted by Sy\(F, G), is an (m + n) x (m + n) matrix constructed from the coefficients of F and 
G, such that 



\ 



Syl(F, G) = f f m g g„ 



\ fo go/ 



(15) 



□ 



Next, we define the "subrcsultant matrix" to derive polynomial subresultants. 



Definition 5 (Subresultant Matrix). Let F and G be defined as in ([14]) . For j < n, the j-th 
subresultant matrix of F and G, denoted by (F, G), is an (m+n— j) x (m + n — 2j) sub-matrix 
ofSy\(F,G) obtained by taking the leftn — j columns of coefficients of F and the leftm — j columns 
of coefficients of G, such that 



N^(F,G) = 



ffn 



V fo 



fo fm go g n 



(16) 



9o J 



□ 



Definition 6 (Subresultant). Let F and G be defined as in (| 14[) . For j < n and k = 0, . . . , j, 

letN { k j) =N { k j) (F,G) be a sub-matrix of N^\F, G) obtained by taking the top m + n — 2j — 1 rows 
and the (m + n — j — k)-th row (note that N ( k j) (F, G) is a square matrix). Then, the polynomial 



Sj{F,G) =det(AfV + 
is called the j-th subresultant of F and G. 



det{N^ } )x° 



(17) 



□ 
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Theorem 1 (Fundamental Theorem of Subresultants |3j). Let F and G be defined as in 
(fT4| . (Pi,...,Pfc) = prs(P, G) be complete PRS, and ((013, $3), . . . , (a.k,(3k)) be its division rule. 
Let m = deg(Pi) and Ci — lc(Pi) for i = 1, . . . , k, and di = rt,; — ni + \ for i = 1, . . ., k — 1. Then, 
we have 

Sj(F,G)=0 forO<j<n k , (18) 
S ni (P,G) ^Pcf'- 1 " 1 ^!^^ 11 I c ^+ di - 1 (-i)(« ! -2-™ l )(« ! -i-« l )| ) ( 19 ) 
Sj(F,G)=0 for m < j < m-i - 1, (20) 

S„ j _ 1 _ 1 (P,G) =^4l*- i n|^J c *-»+«Ji-i(_ 1 )(n,- 9 -n < -i+l)(n,- 1 -n,- 1 +l)| j (2 1) 

for i = 3, . . . , k. □ 

By the above theorem, we can express coefficients of PRS by determinants of matrices whose 
elements are the coefficients of initial polynomials. 



3.2 Recursive Subresultants 

We construct "recursive subresultant matrix" whose determinants represent elements of recursive 
PRS by the coefficients of initial polynomials. 

To make help for readers, we first show an example of recursive subresultant matrix for recursive 
Sturm sequence in Example [1] 

Example 2 (Recursive Subresultant Matrix). We express P(x) and -j^P{x) in Example [T] by 



P(x) = f S x 8 + --- + fox , —F(x) = g r x 7 + ■■■+ g x°. 

dx 



(22) 



Let M {1 ^(F,G) = N {1 ^(F,G), then the matrices M^' 5) (P,G), m£' 5) (F,G) and m'^' 5) (F,G) 
are given as 



M (1 ' 5) (P,G) = 



M, 



(1,5) N 



(h 97 \ 

h fa 96 9i 
h h 95 5e 57 
h h 54 55 56 



M 



h h 53 54 55 

/3 fi 52 53 54 
/2 h 51 52 53 

h h 50 51 52 

/o /1 5o 5i 
Jo 50/ 



M 



'(1,5) 



(F,G) = 



/5/ 4 5/5 553 5c/ 4 5g 5 \ 
4/3 4/4 4 52 4 ff3 4 54 
3/2 3/3 3gi 3g 2 3.93 
2/i 2/a 2 5o 2 3l 2. 92 

V /o /1 5o 5i / 



(23) 



where horizontal lines in matrices divide them into the upper and the lower components. Note 
that the matrix Af'(i> 5 )(p G) is derived from M^' 5 \f,G) by multiplying the 2-th row by 6 — / 
for / = 1, . . . , 5 and deleting the lowest row. Similarly, the (2, 3)-th recursive subresultant matrix 
jVf (2,3) ^ j g cons tructed as 

\ 



M (2 ' 3) (P,G) = 





















0---0 




0---0 



(24) 



□ 
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M (fc ' j) (F,G) = 



/M,r i,,k - ij 



(fc-i,«_i) 



Af, 







M 



(fc-ijfe-i) 



M, 



(fc-ij fc -i) 



Al 



(fe-ijfc-i) 



A/, 



(fc-ijfc-i) 
u 







M 



'{k-Uk-l) 



M, 



(fc-i.j fc -i) 

u 



M 



( l k-\,j k _ l ) 



Fig. 1. Illustration of (F,G). Note that the number of blocks M 

whereas that of M^ k ljfc - 1 - ) is jk-i — j; see Definition [7] for details. 



(fc-ijfc-i) 



is ife-i - j - 1, 



Definition 7 (Recursive Subresultant Matrix). Let _F and G fee defined as in (fT4l) . and Zei 

(P^, ■ ■ ■ , Pj 1 % ■ ■ ■ , Pi* , ■ ■ ■ 7 P/*^) °e complete recursive PRS for F and G as in Definition^ Let 
jo = m and jk = /or fc = 1, . . . , t. Then, for each tuple of numbers (k, j) with k = 1, . . . , t and 
j = jk-i — 2, . . . , 0, define matrix M' fe,J ' (F, G) as follows. 

1. For k = l, let M^(F, G) = N^(F, G). 

2. Fork > 1, let M {k ^(F,G) consist of the upper block and the lower block, defined as follows: 

(a) The upper block is partitioned into (Jk—i —jk — 1) x (jk-l ~ jk ~ 1) blocks with diagonal blocks 
filled with M^ 1 ' jh ~ 1 \ where M^ -1 ^ - ^ is a sub-matrix of Af^ -1 '-"- 1 ) (F, G) obtained by 
deleting the bottom jk-i + 1 rows. 

(b) Let M^- 1 ' 3 "-^ be a sub-matrix of M^ k G) obtained by taking the bottom jk-i + l 

rows, and let M^ 1 be a sub-matrix of ljfc - 1 ) by multiplying the (jk-i + I — t)- 

th rows by r for r = jk-%, ■ ■■,!> then by deleting the bottom row. Then, the lower block 

consists of jk-i — j — 1 blocks of l,3h ~ 1 ' such that the leftmost block is placed at the 

top row of the container block and the right-side block is placed down by 1 row from the 

1 (k i j t _i ) 

left-side block, then followed by jk-i — j blocks of placed by the same manner 

as M^- 1Jk - l} . 

As a result, (F, G) becomes as shown in Fig.{J\ Then, M^ k,j \F,G) is called the (k,j)-th 

recursive subresultant matrix of F and G. □ 

Proposition 1. Fork = 1, . . . ,t and j < jk-i—1, the numbers of rows and columns of M^ k '^ (F, G), 
the (k,j)-th recursive subresultant matrix of F and G are 

(m + n - 2j x ) | J[ (2j,_ x - 2j t - 1) j {2j k -i - 2j - 1) + j (25) 

and 
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respectively. 

Proof. By induction on k. For k = 1, by definition of the subresultant matrix, we see that 
MV-'i>{F, G) is a (to + n — j) x (m + n — 2j) matrix. Let us assume that equations ([2"5)) and ([25]) 
are valid for 1, . . . , k — 1. Then, we calculate the numbers of the rows and columns of M<- k ^ (F, G) 
as follows. 

1. The numbers of rows of ljfc ^^ an d fe_1 '' 7fc ~ 1 ' are equal to jk-i + 1 and j k -i, respec- 
tively, hence the number of rows a block which consists of M^ fc and M^* 1,Jh - 1 ^ i n 
M^(F,G) equals 

2j*-i-j-l. (27) 

On the other hand, the number of rows of 1j 'js-i) is equal to (m + n — 2ji){nf=2 1 (^-i — 
2ji — 1)} — 1) hence the number of rows of diagonal blocks in M^ k \F, G) is equal to 

(m + n - 2jr) j [] (2j,_i " 2j; - 1) - l| (2j fe _x - 2j - 1). (28) 

By adding formulas ([27]) and ([28)). we obtain |25)) . 

2. The number of columns of M^ -1 '- 7 '* -1 ^^ G) is equal to (m + n - 2ji){nfr 2 1 (2j ; _i - 2j ; - 1)}, 
hence the number of columns of M^ k '^(F, G) is equal to ([26]) . 

This proves the proposition. □ 

Now, we define recursive subresultants by recursive subresultant matrices. 

Definition 8 (Recursive Subresultant). Let F and G be defined as in (lbl|) . and let (P x (1) , . . . , 
P^\ ...,P^\..., P// 3 ) 6e complete recursive PRS for F and G as in Definition^ Let jo = m and 

jk = for k = 1, . . . , t. For j = j k _ x - 2, . . . , and r = j, . . . , 0, let M ( T k ' j) = M?' j) (F, G) be 
a sub-matrix of the (k,j)-th recursive subresultant matrix M^ k '^ (F,G) obtained by taking the top 
(TO + n-2j 1 ){nfr 2 1 (2j;-i-2j i -l)}(2j fc _ 1 -2j-l)-l rows and the ((m + n-2j 1 ){l\^2 - 
2ji — i)}(2jk—i — 2j—i)+j — T)-th row (note that is a square matrix). Then, the polynomial 

S Kj (F, G) = det(Mj fcj V + ' ' ' + det(M (fcJ) )a; (29) 
is called the (k, j)-th recursive subresultant of F and G. □ 
Finally, we derive the relation between recursive subresultants and coefficients in recursive PRS. 

Lemma 1. Let F and G be defined as in (|14[) . and let (Pi , ■ . ■ , P[ l \ • ■ ■ , p[ 1 \ ■ ■ ■ , Pj^) be com- 
plete recursive PRS for F and G as in Definition^ Let c^ — lc(P^), — deg(P^), jo = rn 

and jk — Tii f or ^ = 1? • • • ! ^ an d * = 1) • • • j lk> an d let d^ = — r^+i f or k — 1, ■ ■ ■ , i and 
i = 1, . . . , Ik — 1. Furthermore, for k = 1, . . . , t — 1 and j = jk-i — 2, . . . , 0, define 

u ktj =(m + n- 2ji) { JJ(2jz-i - 2j; -1)1 (2j fc _i - 2j - 1), 



i=2 




(30) 



s fc = (O II Hit (c^x)^-^-^^!)^ 



and let Uk = Uk,j k ■ For k = 2, . . . , t and j = jk-i — 2, . . . , 0, define 

b k<j = 2 Jfc „x - 2j - 1, r kJ = ( _i ) K- 1 -i)(i+2 + ...+(6 fe ,,-i)) ! (31) 

and Zet = b ky j h and r k = rfcj fc . Then, for the (k,j)-th recursive subresultant of F and G with 
k = 1, . . . ,t and j = jk-i — 2, . . . , 7 we have 

S k>j (F, G) = R ktj ■ SjiP^lpM), (32) 

where R ltj = 1 and R k>j = ((••• ({B^ ■ r 2 B 2 ) b - ■ r 3 B 3 ) b * ■ ■ ■ ) b ^ ■ r k -iB k -i) b ^ ■ r kJ for k > 1. 
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Proof. For a univariate polynomial P(x) = a n x n + • ■ ■ + oqx with aj G R for j = 0, . . . , n, let us 
denote the coefficient vector for P{x) by p = t (a n , . . . , a ). 

We prove the lemma by induction on k. For k = 1, it follows immediately from the Fundamental 
Theorem of subresultants (Theorem [T). Let us assume that (|3"2|) is valid for 1, . . . , k— 1. Then, we 
have 

S^^^PG) = R* . • S^fPf-Vjj*- 1 '), (33) 

hence we have 



det(M^-^-i)) = Bfc-ij,., •det(^-)(P 1 (fc - 1) ,P 2 (fc - 1) )), 



(34) 



for r = jfc-i, ... ,0. Therefore, there exists a matrix Wu-i such that det(Wfc-i) = Rk-i,j k _ 1 and 
that we can transform M^ -1 ^- 1 ' (P, G) to 



JV«*-i)(Pi (fc_1, ,P 2 tfc " 1J ) / ' 



(35) 



by eliminations on columns. Furthermore, by eliminations and exchanges on columns 
in the block N^k-i) (p( k x \p^ k as shown in Brown and Traub ([3]), we can transform 
M&-ijh-i)(F,G) to 



iV, 



(Jfc-l) 

C7 



V 



Pi 



(36) 



such that iV, 



(ifc-i) 



is a lower triangular matrix with all diagonal elements equal to 1 



and that det (Mr 



>(k-l,jfc-i) 



(fc-i j'fc-i) 



(P,G)) = P fe _i ■ det(M)' 



(fe-ijfc-i) 



(P,G)), where Aff 



(fc-i jfc-i) 



(PG) 



(P, G) are sub-matrices of M^" 1 (P, G) and M^' 1 ^"- 1 \F, G)), respectively, 



and M) 

obtained by the same manner as we have obtained M {k (F, G). Therefore, we have 

det(M^- ljfc - 1 )(P,G)) = Pfe_i • det(A 7 4 fc - ljfc - l) (P,G)). 



(37) 



Similarly, let M'( fc_1 ' ife - 1 )(P, G) 



(fc-ijfc_i) N 
c; 

'(k-l, 3k -l) 



. Then, by the same transformations in the 



above, we can transform M ^ k 1 >-"=- I )(P ) G) to 







o 














{ * 







(38) 



satisfying 



det(M T ' 



/(fe-i,ife-i) 



(P,G)) = P fc _! -det(M T ! 



'(fe-ij fe _i) 



(PG)), 



(39) 

where M^^-'^P, G) and M T {k ~ 1 ' jk - l) (F,G)) are sub-matrices of m'^ 1 -""- 1 ) (P, G) and 
M (.k-ijk-i) (jr, G)), respectively, obtained by taking the top Uk-i — 1 rows and the (tik-i + 
— r)-th row for r = jk-i . .., 1. Therefore, for j < jk-i — 1, we can transform M^ k ^\F, G) 
to M( k 'i\F,G) as shown in Fig. [2] by eliminations and exchanges on columns in each column 
block, and let M? ,j) (P, G) be sub-matrix of M {k ^ (P, G) obtained by the same manner as we 
have obtained M^ fe,i) (P,G). Then, we have 



det(M^(PG)) - (B k -i) bk - ■ det(M T k ^(F,G)), 



(40) 



from (37) and (39) and since there exist & fc)j blocks of M^ -1 '- 7 '*- 1 )^, anc i Af'^ -1 '- 7 '*-^^, G) in 
M( k 'ti(F, G) with each of which divided into the upper and the lower block. 

Furthermore, by exchanges on columns, we can transform M {k ^ (P, G) to M {k ^ (P, G) as shown 
in Fig. [3j and let M^'^ (P G) be sub-matrix of Af( fe > J ) (P, G) obtained by the same manner as we 
have obtained M^ k,j) (F,G). Then, we have 
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M (k ' j) (F,G) = 





6 

Jv r/ 





















o 









»'fe-l 


Jv rr 













0-0 

































\ 




o 











o 











o 







_(*) 


0-0 










_(*) 

P2 












(fc) 

P2 / 



Fig. 2. Illustration of M^\F, G). See Lemma H]for details. 



det(M^ fc ^ (F, G)) = r fc j • det(M<* ^ (F, G)), (41) 



because the (ukj — (I — l)iifc_i)-th column in M^ k '^(F, G) was moved to the (v,k,j — (I — l))-th 
column in M^ k, ^(F, G) for 1 = 1,..., bkj- Furthermore, we have 



det(M( fc ^(F,G)) = (R h ,., : lh. L ) 6) • N^{p[ k \p^ k) ), (42) 



because the lower-right block of p[ k) and pf> in M^\F, G) is equal to N^(Pf\ P 2 (fe) ). 
Finally, from flD]). (|H|) and (gSJ), we have 



det(M^(F,G)) = r kij ■ (R^j^B^)"^ ■ det(JV«>(if \ P<" >)) 
= i? fej -det(7V«(P 1 (fe) ,P 2 (fc) )). 



Therefore, by the definitions of recursive subresultant, we obtain (|32p . This proves the lemma. □ 
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M {k ' 3) {F,G) = 



-111 



111 



hi 



111 



(k) 



(fc) 

*>2 



p 2 



„(*) 



Fig. 3. Illustration of M^'fi (F, G). Note that the lower-right block which consists of pj" and P2 
is equal to N^(P^ k) , P 2 (fe) ), and the number of blocks Wfc_i and iV^ _1) is b fc j = 2j fc _i - 2j - 1: 
see Lemma [1] for details. 



Theorem 2. M^zt/i £/ie same conditions as in Lemma{^ and for k = 1, . . . , t and i = 3,4, . . . , Ik, 
we have 



S ktj (F,G)=0 forO<j <ng\ 
S k ^ ) (F,G) = P^(4 k) ) d ^- 1 R 



fe,-. 



n 

;=3 



^f 1 ' 1 "* ( c W i )^-U^ , 1 )(_i)K < - ) 2 -™^ , )K ( - , 1 -«< fc) ) y 



S k<j (F, G) = for n[ k) < j < nf\ - 1, 

s fc ,„«- 1 (^G)=^( C fj 1 ) i -e\ Ji 



n 

;=3 



(fe) \ n i-i 



n, , —n) 



( c wj(4: j 2+4- ; i)(_i)(< J 2-^- ; 1 +i)K c - ) x-»& ) 



a 



i+l) 



(44) 

(45) 
(46) 

(47) 



Proof. By substituting S^Pf , P 2 W ) in (32]) by ip ]H [2"I |) . we obtain pi |) -((17 |l . respectively. □ 
We show an example of the proof of Lemma [1] for recursive subresultant matrix in Example 
Example 3. Let us express P^ in Example [T] by 



(48) 
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with rij = deg(P^ > ). By eliminations and exchanges of columns as shown in Brown and Traub ([3]), 



,(fc) 



we can transform M^ b \F,G) = 
and M ^^(F, G), respectively, as 



SandM'(^)(F,G) 



M, 



(1,5)' 



, M 



T (XW 



to M^(F,G) 



M (1,5) (F,G) 



N, 



m' (1 ' s )(F,G) 



(5) 



4 5) 





* 


pf 



a (1) 

"2,6 

"2,5 
(1) 
2,4 



"2,3 

"2,2 
(1) 
2.1 
(1) 
2.0 



1 

f 2,6 
l 2,5 



a 



J 17 

"2,4 

a (1) 

"2,3 

a (1) 

"2,2 

a (1) 

"2,1 

(1) 



(1) 
3,5 



1 



(I 



J 17 

a 3,4 

a (1) 

"3,3 

a (1) 

"3,2 

a (1) 

"3,1 

a (1) 



2,0 "3,0 



"3,4 

a (1) 

"3,3 

a (1) 

"3,2 

a (1) 

"3,1 
(1) 
3,0 



a 



"4,4 

"4,3 

a (1) 
"4,2 

a (1) 

"4,1 
(1) 



1 

(1) 
2.6 
(1) 
2,5 
(1) 
2 4 

5a 2 3 
4a« 
3«g 
2a (1) 

z "2,0 



1 

a (1) 
"2,6 

(i) 



2.5 



a 



(1) 
3.5 



4,0 ' 



1 



°"3,3 

2a 3 2 
(i) 



2.4 
3a (1) 

°"2 2 

(i 



3,4 
3a (1) 

oa 3,2 

(i) 



2.0 



a 



3,0 



3,1 



4a 



°"4,5 
(1) 
4.4 

3a (1) 

4!ly 



(49) 



where 4 j = af-j /a^l- Furthermore, we have 



det(M^' 5 )(F, G)) = B x • det(M^< 5 )(F, G)) for r = 5, . . . , 0, 
det{M^ lfi \F, G)) = B x ■ det(M T ^' 5 \F, G)) for r = 5, . . . , 1, 



(50) 



with 



Bx - -(a, 7 ) (a 3 6 ) , 



(51) 



where A#' 5) (F, G) and m' t {1,5) (F, G)) are sub-matrices of M^^(F, G) and M'^^(F, G), respec- 
tively, obtained by taking the top 4 rows and the (10 — r)-th row. Therefore, by eliminations and 
exchanges on columns, we can transform M^ 2,3 \F, G) in (f2"4")l to M^ 2 ' 3 \F,G) as 



M (2 ' 3) (F,G) = 



(2) 
Pi 



N, 



(2) 
P 2 



0---0 



ivjn 



0---0 



(2) 

p\ 



(52) 



374 Akira Terui 



satisfying det(M^ 2 ' 3) (P, G)) = (Pi) 3 • det(M^ 2 ' 3) (P, G)). Furthermore, by exchan ges on columns, 
we can transform M^ 2 ' 3 \F, G) to Af( 2 ' 3 >(P, G) as 



M (2 ' 3) (P,G) = 



\ 



0---0 



N, 



0---0 



(2; 
Pi 



(2) 

p\ 



(2) 
P\ 



(53) 



N, 



u 



N, 



satisfying det(M^ 2 ' 3) (F, G)) = r 2 , 3 • det(M^ 2 ' 3) (F, G)) = r 2 , 3 • det(ivf ) (F 1 (2) , P 2 (2) )). Therefore, we 
have 

det(M^ 3 )(F, GO) = (i?i) 3 r 2 , 3 ■ det(AK 3 )(P 1 (2) , P 2 (2) )) = P 2 , 3 • det^CP^, P 2 (2) )), (54) 
for r = 3, . . . , 0, and we have 

S 2 , 3 (P,G) = P. 3 • SsCPf.Pf) = {(ag) 2 (a«) 2 } 3 (45) 2 ■ Pf . (55) 

□ 

4 Conclusion and Motivation 

In this paper, we have defined recursive PRS as well as recursive subresultants, and proved a similar 
theorem as the fundamental theorem of subresultant. 

The concept of recursive subresultant is inspired, in approximate algebraic computation, by 
representing coefficients in recursive PRS by the coefficients of initial polynomials. For example, 
consider calculating recursive Sturm sequence of a polynomial with floating-point number coeffi- 
cients by floating-point arithmetic. In the case the initial polynomial has multiple or close zeros, 
there may exist a polynomial in the sequence such that it is difficult to decide whether the poly- 
nomial becomes zero or not. Also, zero recognition of very small leading coefficient is another 
important problem because it plays crucial role in calculating the number of real zeros. 

For the problem of zero recognition of very small leading coefficients, the present author and 
Sasaki ([5]) have proposed a criterion for calculating the number of real zeros correctly by Sturm's 
method: if the Sturm sequence satisfy certain condition on Sylvester matrix, then we can neglect 
the small leading coefficient which makes computation of the Sturm sequence more stable. We 
expect that the recursive subresultant (matrix) will be useful for zero recognition of a polynomial 
in recursive Sturm sequence, by representing its coefficients by the coefficients of initial polynomials 
then analyzing it by numerical methods; this is the problem on which we are working now. 
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